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ABSTRACT 

We compare the B-band luminosity function of virialized halos with the mass function predicted by 
the Press-Schechter theory in cold dark matter cosmogonies. We find that all cosmological models fail to 
match our results if a constant mass-to-light ratio is assumed. In order for these models to match the 
faint end of the luminosity function, a mass-to-light ratio decreasing with luminosity as £~°- 5±0 - 06 is 
required. For a ACDM model, the mass-to-light function has a minimum of ~ in solar units in 

the i?-band, corresponding to ~ 25% of the baryons in the form of stars, and this minimum occurs close 
to the luminosity of an L* galaxy. At the high-mass end, the ACDM model requires a mass-to-light 
ratio increasing with luminosity as r,+°- 5±0 - 26 . This scaling behavior of the mass-to-light ratio appears 
to be in qualitative agreement with the predictions of semi-analytical models of galaxy formation. In 
contrast, for the rCDM model, a constant mass-to-light ratio suffices to match the high-mass end. 

We also derive the halo occupation number, i.e. the number of galaxies brighter than L* al hosted in 
a virialized system. We find that the halo occupation number scales non-linearly with the total mass of 
the system, N gai {> L* gal ) oc m o.55±o.026 for the ACDM model. 

We find a break in the power-law slope of the X-ray-to-optical luminosity relation, independent of the 
cosmological model. This break occurs at a scale corresponding to poor groups. In the ACDM model, 
the poor-group mass is also the scale at which the mass-to-light ratio of virialized systems begins to 
increase. This correspondence suggests a physical link between star formation and the X-ray properties 
of halos, possibly due to preheating by supernovae or to efficient cooling of low-entropy gas into galaxies. 

Subject headings: cosmology: large-scale structure of the universe — cosmology: dark matter — 

galaxies: clusters: general — galaxies: halos — galaxies: luminosity function, mass 
function — X-rays: general 



1. introduction 

One of the major problems of cosmology is to determine 
the connection between mass and light in the universe over 
different scales. The abundance by mass of virialized dark 
matter halos is a fundamental prediction of cosmological 
models (Press & Schechter 1974). Observationally, halos 
must be identified with virialized galaxy "systems" , rang- 
ing in size from single dwarf galaxies to rich clusters of 
galaxies. 

The masses, and hence the mass multiplicity functions, 
derived from cataloged systems are not robust, particu- 
larly for small systems such as poor groups which contain 
only a few galaxies. In contrast, halo light is a funda- 
mental observable obeying a robustly determined distribu- 
tion function, the luminosity function of virialized systems 
(VSLF). 

In a previous paper (Marinoni et al. 2001b, Paper V), 
we measured the VSLF in the nearby Universe using an 
objective catalog of virialized systems. In this paper we 
show how the VSLF can be used to investigate scaling 
relations between mass and light. 

The direct approach to obtaining scaling relations be- 
tween, for example, mass and light is to regress one param- 
eter on the other. This direct approach has several weak- 
nesses, however. As noted above, for low mass systems 



the measured masses can have large errors. Second, the 
catalogs may be incomplete in either mass or light. Fur- 
thermore, it is usually necessary to patch together results 
from several independent catalogs, each of which probes 
different mass scales (e.g., galaxies, groups and clusters). 

In this paper we take a different approach. We derive 
the mass functions from theory and compare these with 
the observed VSLF, solving for the unknown mass-to- 
light ratio (T) as a function of the total mass or lumi- 
nosity of the halo. The price paid in this approach is that 
the mass function depends on the assumed cosmological 
model. The advantage over the conventional approach is 
that we can simultaneously and seamlessly explore a wide 
dynamic range in mass and luminosity, free of systematics. 

The T function, which measures the efficiency with 
which the universe transforms matter into light, can thus 
be used not only as a traditional estimator of the density 
parameter f2 m but also as a diagnostic tool for different 
models of galaxy formation. 

A closely related statistic, useful in modeling galaxy for- 
mation and clustering, is the halo occupation number, i.e. 
the number of galaxies above some luminosity threshold 
(chosen to be the absolute magnitude of an L* galaxy) 
that populate a virialized halo as a function of its total 
mass or luminosity. 
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Within the same framework, we can also investigate the 
emission properties of halos in different wavelength do- 
mains of the electromagnetic spectrum such as in the op- 
tical and X-ray bands. The sample of systems for which 
both the optical and X-ray luminosities are well measured 
is still limited, and does not allow a reliable determination 
of the scaling of the X-ray to optical light ratio. However, 
by comparing statistically their luminosity distributions, 
we can constrain the functional behavior of this ratio. 

This paper is the sixth in a series (Marinoni et al. 1998, 
Paper I; Marinoni et al. 1999, Paper II; Giuricin et al. 2000, 
Paper III; Giuricin et al. 2001, Paper IV; Marinoni et al. 
2001b, Paper V) in which we investigate the properties of 
the large-scale structures as traced by the NOG sample. 
The outline of our paper is as follows: in §2 we review and 
summarize the determination of the observed luminosity 
function of virializcd halos. In §3 we compare our results 
with the predictions of rCDM and ACDM models, as well 
as the predictions from semi-analytic models. In §4, we 
measure the halo occupation number. In §5 we examine 
the relationship between the optical and X-ray properties 
in virialized systems. In §6 we discuss in detail the result- 
ing scaling relations between the properties of groups and 
clusters of galaxies. Results are summarized in §7. 

Throughout this paper, the Hubble constant is taken to 
be 75 /175 km s _1 Mpc -1 and the recession velocities czlg 
are evaluated in the Local Group rest frame. 

2. THE LUMINOSITY FUNCTION OF VIRIALIZED HALOS 

In paper V, we determined the luminosity function of 
virialized systems using groups extracted from the NOG 
galaxy catalog. The NOG sample (Marinoni 2001a) is 
a statistically controlled, distance-limited (czlg <6000 
km s _1 ) and magnitude-limited (B < 14) complete sam- 
ple of more than 7000 optical galaxies. The sample cov- 
ers 2/3 (8.27 sr) of the sky (|6| > 20°), a volume of 
1.41 x 10 6 /iy 5 Mpc~ 3 and has a redshift completeness of 
98%. 

From the NOG, three different "virialized system" sam- 
ples were determined using groups selected by three dif- 
ferent objective group- finding algorithms (Paper III). In 
paper V, we showed that the VSLF derived from these 
different samples are consistent and hence that the VSLF 
is robust. 

In order to recover accurately the absolute luminosity of 
a given halo, we correct for both the integrated luminosity 
of galaxies below the magnitude limit of the sample and for 
the completeness of subsamples of the halo catalog. More- 
over, models to correct for the large-scale motions (Paper 
I) have been applied in order to improve the determination 
of the halo distances. The key steps of the reconstruction 
procedure are shown in Fig. 1, focusing on several virial- 
ized systems in the Virgo region which are embedded in 
lower-density filamentary large-scale structure. 

In Paper V, we found that the £>-band luminosity func- 
tion of the whole hierarchy of gravitationally bound sys- 
tems, from single galaxies to rich clusters, was insensi- 
tive to the choice of the group-finding algorithms with 
which halos are selected, and was well described over 
the absolute-magnitude range —24.5 < M s + 5 log ^75 < 
— 18.5 by a Schechter function (Schechter 1976) with a s = 
-1.4 ± 0.03, M* - 51og/i 75 = -23.1 ± 0.06 and cj>* s = 



4.8 x 10 4 hy 5 Mpc 3 or by a double power law: (f> p i(L s ) oc 

£ -1.46±0.07 for Lg < Lpl and 0( is ) K L -2.35f0.16 for 

L s > Lpi with Lpi — 8.5 x 10 10 /i 75 2 L Q , corresponding to 
M s — 51og/i75 = —21.85. The characteristic luminosity 
of virialized systems, L p \, is ~ 3 times brighter than that 
(I/* al ) of the luminosity function of NOG galaxies. 

3. THE MASS-TO-LIGHT RATIO OF VIRIALIZED SYSTEMS 

3.1. Introduction 

In order to compare observational data with theoreti- 
cal predictions, the traditional approach has been to esti- 
mate the halo mass of groups and clusters directly from 
the data, using, for example, velocity dispersions of galaxy 
members (Carlberg et al. 1996; Girardi et al. 1998), X-ray 
gas temperatures (David, Jones, & Forman 1996; Lewis 
et al. 1999), and gravitational lensing (Smail et al. 1997, 
Allen 1998). 

The application of these methods is reasonably straight- 
forward for rich clusters of galaxies but is more problem- 
atic for poor groups, due to their low X-ray surface bright- 
ness and the poor sampling of optical dynamical mass esti- 
mators (see, e.g., Zabludoff & Mulchaey 1998, Mahdavi et 
al. 1999). Consider, for example, a system composed of a 
few elements selected from a magnitude-limited survey. In 
this case one will miss faint members that do not satisfy 
selection criteria. The contribution of faint members to 
the total luminosity of the system is negligible small, and 
can easily be corrected (see Paper V). In contrast, these 
faint galaxies are important dynamical probes of the grav- 
itational potential well, and their omission can translate 
into a serious bias in the estimated mass, particularly if 
galaxies of different luminosities are clustered in different 
ways (Park et al. 1994; Giuricin et al. 2001). Moreover 
one can construct a completeness function and than derive 
a differential distribution function for the light emitted by 
virialized systems, while the construction of a similar com- 
pleteness function involving dynamical estimators (veloc- 
ity dispersion for example) is very problematic (Borgani 
et al. 1997). Thus, the luminosity function of virializcd 
systems is more robust on group scales than is the mass 
function derived using projected velocity dispersions. 

In this paper, following Cavaliere, Colafrancesco & 
Scaramella (1991) and Moore, Frcnk, & White (1993, 
MFW), we explore an alternative approach to determin- 
ing the mass-to-light ratio over a range of scales including 
both poor and rich systems. Specifically, we compare the 
robust observationally-determined luminosity function of 
systems to theoretically-determined mass functions. 

3.2. Theory 

The differential mass function of virialized systems n(m) 
can be calculated, for a given cosmological model, via the 
Press-Schcchter formalism (Press & Schechter 1974, here- 
after PS; Schaeffer & Silk 1985; Bower 1991; Bond et al. 
1991; Cavaliere, Colafrancesco, & Scaramella 1991; Lacey 
& Cole 1993; Monaco 1998), or more accurately, through 
N-body simulations (Efstathiou & Rees 1988; Lacey & 
Cole 1994; Governato et al. 1999; Jenkins et al. 2001). The 
comparison with N-body simulations reveals that the PS 
mass function provides a reasonably accurate description 
of the abundance of virialized halos on group and cluster 
scales although it tends to overpredict 
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Fig. 1. — Upper left: Density (5=1.5) distribution of NOG galax- 
ies smoothed using a Gaussian window function with smoothing 
length R s = 200 km s _1 . Upper right: redshift distribution of NOG 
galaxies in a cone diagram centered on the Virgo region (smoothed 
with _R S = 100 km s _1 ). Regions with <5 > are in grey scale while 
contours are spaced by A<5 = 2 Lower left: the same region af- 
ter clustering reconstruction with the hierarchical method. (Lower 
right: the same region after having applied the peculiar velocity 
field model fitted using Mark III data. 
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(albeit by a small amount) the abundance of low-mass 
halos and underpredict that of high-mass halos (e.g., Gross 
et al. 1998; Sheth & Tormcn 1998; Governato et al. 1999). 
Given the uncertainties in the luminosity function of sys- 
tems to which we will be comparing the mass function, the 
PS formula is of sufficient accuracy. 

The PS analytical expression for the present-day comov- 
ing number density of dark halos of mass m in the interval 
dm is 



n(m)dm 



pS c 



7r a(m) 2 m 



da(m) 
dm 



ex P I ~ n _/ C _N9 ) ( ) dm 

(1) 

where p is the present mean density of the universe, a(m) 
is the present-day rms fluctuations of the linear density 
field after smoothing with a spherical top-hat filter con- 
taining a mean mass m (and it is specified by the power 
spectrum, P(k), of the density fluctuations), and S c is a 
density threshold usually taken to be the extrapolated lin- 
ear overdensity of a spherical perturbation at the time it 
collapses. The parameter S c is only weakly-dependent on 
matter density parameter f2 m and the cosmological con- 
stant ft a (varying by less than 0.02 for flat models with 
n m > 0.1, see Fig 1. of Eke, Cole, & Frenk 1996), and 
so has been set to its Einsten-de Sitter cosmology value 
1.686. 

The dimcnsionless power spectrum (A(fc) = 
t^(27r 2 ) _1 fc 3 |(5fc| 2 in the linear regime is given by the ana- 
lytic approximation of Bond & Efstathiou (1984) 



A(fc) 



Ak A 



1 + [aq + (6g)3/ 2 + (eg) 2 ]" 



2/v 



(2) 



where q 
Mpc, c 



fc/r, a = 6.4/i 10 o Mpc, 



b = 3fc7/oo 



= 1.7ft 



100 



Mpc, v = 1.13, and h w $ = 
i?o/(100km s- 1 Mpc -1 ). The parameter T = n m h 10 o (if 
we neglect the effect of baryons) determines the shape, and 
the normalization is determined e.g. by Cosmic Microwave 
Background fluctuations or by fixing the rms fluctuations 
on a scale of 8ftf Q Mpc, er 8 . 



3.3. Results 

Figure 2 shows the PS-predicted luminosity functions 
obtained assuming constant mass-to-light ratio for four 
different cosmological models: Standard CDM (O m = 1, 
n A = 0, H = 50 km s" 1 Mpc" 1 , er 8 = 0.51), ACDM 
(fi m = 0.3, n A = 0.7, H = 70 km s" 1 Mpc -1 , er 8 = 0.9), 
rCDM (O m = 1, H = 50 km s -1 Mpc" 1 , cr 8 = 0.51, T = 
0.21) and Open CDM (fi m = 0.3, n A = 0, H = 70km s -1 
Mpc -1 , erg = 0.85). These models are described in detail 
by Jenkins et al. (1998) and are normalized adopting the 
values for the fluctuations in an 8 ft^~ 00 Mpc sphere, cr 8 , 
prescribed by Eke, Cole, & Frenk (1996) from their anal- 
ysis of the local cluster X-ray temperature function. 

The NOG VSLF and the corresponding MFW results 
are shown in Figure 2. The latter is represented by a 
hatched region, which shows how the VSLF varies if the 
overdensity contrast used in identifying groups is varied 
from 50 to 300. 
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Fig. 2. — The NOG and CfA system LFs are compared with the 
luminosity distribution predicted using the Press-Schechter func- 
tion and a constant mass-to— light ratio. Four different cosmological 
models (the standard CDM, ACDM tCDM and open CDM) are 
computed with cosmological parameters as given in the figure. For 
each model we show the constant mass— to-light ratio adopted. The 
hatched region indicates how the LF varies if the overdensity crite- 
rion — used in identifying groups is varied from 50 to 300. 



The friends-of-friends grouping algorithm used to define 
the groups and hence the VSLF requires that the local den- 
sity contrast at the edge of a group be greater than some 
limiting threshold ^ (see the discussion in §2). If we as- 
sume a spherical halo with a singular isothermal density 
profile n(r) oc r~ 2 , this local density threshold corresponds 
to a mean overdensity {^-) — 3^ where the average is 
performed over the typical separation of group members 
(virialization region). Thus, if galaxy biasing is negligible, 
for our grouping algorithms, this implies a value which is 
close to the mean overdensity of nearly 180 predicted for a 
top-hat spherical collapse in the case of Cl m = 1 (see Eke, 
Cole, & Frenk 1996 for values in different cosmologies). 

Cole & Lacey (1996) found that the virialized regions 
of N-body halos are reliably identified using the percola- 
tion method and a local density threshold parameter of 
nearly 60. Using this limiting density contrast, they are 
able to reproduce the bulk properties of virialized struc- 
tures (in particular the virial mass) and concluded that the 
closeness to global virial equilibrium of the identified N- 
body halos depends rather weakly on the adopted limiting 
threshold. Moreover, as noted by Diaferio et al. (1999), 
the average luminosities of groups objectively identified 
from a redshift survey with typical friends-of-friends pa- 
rameters are in agreement with those of groups extracted 
from real-space simulations. Thus, we can use our VSLF 
for a meaningful comparison with the PS predictions. 

Figure 2 indicates that, in agreement with previous re- 
sults, all models fail to describe the observed faint-end 
slope. 
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Fig. 3. — The NOG and CfA system LFs arc compared with the lu- 
minosity distribution predicted using the Press-Schechter function. 
The LF in the tCDM and ACDM scenarios are computed for con- 
stant mass-to-light ratios (T = 33O/175T0 and T = 96/175X9 re- 
spectively) and also (heavy lines) using the piecewise mass-to-light 
ratios given in the figure. The hatched region indicates how the LF 
varies if the overdensity criterion £2 used in identifying groups is 
varied from 50 to 300. 

More interestingly, several models (ACDM and OCDM) 
fail to reproduce the bright end of the LF, if a constant 
mass-to-light ratio is assumed. 

We can turn this problem around, and solve for the 
mass-to-light function, T(m), required to fit the ob- 
served VSLF. We assume that there is a monotonic re- 
lation between the mass and light of virialized systems, 
i.e. L s = ./(to) with no scatter. 

It is possible to obtain an exact numerical solution by 
noting that the T function is defined implicitly by the 
requirement that the number density of systems above 
with masses greater than mass m must be the same as 
the number density of systems with luminosities greater 
than L s = f(m), 

N s {> L s = f(m)]=N s (>m) (3) 

where iV s (> L s ) is given by the appropriate integral over 
the VSLF 

N s (> L s ) = / cf>(L s )dL s (4) 
Jl s 

and N s (> to) is the corresponding integral over the PS 
mass function. 

The LF integral can be obtained analytically, but the 
PS mass function integral must be solved numerically. The 
mass-to-light ratio is then given by 



M s -5 log h 75 

-16 -18 -20 -22 -24 



4 






















tCDM exact 


3 








2 




ACDM 


tCDM 




_ T = 8.6 10" (_L)- Q 


5 h T& T e if M e >-21. 


T=4.8 10' (JO"" h TS T u if M„>-21.85 _ 


1 


_ T = 96 h 75 T Q 


if -23.5<M s <-21 


T=330 h T5 T G if M^-21.85 _ 




_ T=2.7 lO" 4 (10° 


5 h J& T if M 3 <-23.5 





8 9 10 11 12 13 



Fig. 4. — Mass— to-light ratio as a function of luminosity for the 
ACDM and rCDM models. The thin solid and dotted lines indicate 
the mass-to— light ratios for the piecewise models shown in Fig. 3 for 
the ACDM and rCDM models respectively, whereas the thick lines 
are the exact solutions obtained by solving eq. 3. The dashed line 
indicates the critical value of the mass-to— light ratio. The hatched 
region represents extrapolation to regions not covered by our data. 

If we assume that mass and light are simply related via 
a power law expression T oc LJ, the halo luminosity func- 
tion assumes the simple analytical form 

<t>{M s ) oc mn(m)(l + 7) (6) 

The predicted VSLFs obtained fitting "piecewise" T 
functions are shown in Figure 3 for the tCDM and ACDM 
models, which are used here as our reference models. 

In Table 2, we give the masses assigned to luminous 
systems for the numerical solution of the T function. In 
Figure 4, the piecewise and numerical solutions for T are 
plotted as a function of luminosity. The result is that, for 
both models, T must go approximately as £,T - 5±0 06 on 
small scales. Interestingly, this behavior agrees with the 
direct determinations of T derived for spiral galaxies from 
the analysis of their rotation curves and for ellipticals from 
a variety of tracers of galactic gravitational field (see the 
review by Salucci & Persic 1997). 

The T function reaches a minimum at a luminosity of 
approximately 10 10 5 L Q , which corresponds roughly 
to L* al . Then it remains quite flat over ~ 1 dcx in 
luminosity. The behavior at the bright end is model- 
dependent: for the tCDM model, T must remain con- 
stant (to oc Lg whereas for the ACDM model it must 
increase with luminosity as £°- 5±0 - 26 from the scale of 
galaxy groups (to ~ I0 13 to q h^ 5 ) to that of rich clusters 

( ~ IO 15 TO h^). 
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Fig. 5. — Upper: the NOG VSLF (points with ±1<t error bars) 
is shown in comparison with the semi-analytic prediction of Ben- 
son ct al. (2000). Results from their rCDM and ACDM models are 
plotted with short-dashed lines and solid lines respectively. Lower: 
mass— to-light ratios from Fig. 4 (smooth curves) are compared to 
those derived by Benson et al. (curves with error bars). 

Recall that our mass-to-light ratios are based on the PS 
mass functions. If, as mentioned above, the PS prescrip- 
tion underestimates the abundance of high-mass objects, 
the true scaling at the high-luminosity end is likely to be 
somewhat steeper than found above. This would also be 
true if we used a lower value S c = 1.3 in the PS mass func- 
tion, as might be expected for a non-spherical collapse for 
bound structures. 

3.4. Discussion 

It is interesting to compare the T functions derived in 
the previous subsection with values found via direct mass 
estimates. As discussed above, these direct methods are 
problematic for low mass systems such as single galaxies 
(because the full extent of the dark matter halo is not 
probed) or poor groups (because of the small numbers 
of dynamical test particles). The T values found in the 
literature are consistent with the range <~ IOO/175Y0 to 
- 3OO/i 75 T spanned by the ACDM and tCDM models 
considered here. 

For high masses, we can compare our results with the 
scaling behavior of T obtained via direct determinations of 
the virial masses of galaxy clusters. Schacffcr ct al. (1993) 
found that Ty cx L y 3 . Girardi et al. (2000) analyzed 105 
rich clusters and found T^ , oc L° B 2 ~ 3 . 

They also noted that a variation of T with scale could 
not be explained by a higher fraction of spirals in poorer 
clusters, thus suggesting that a similar result would also 
be found by using R-band galaxy magnitudes. Carlberg 
et al. (1996) found that the number of bright galaxies per 
unit mass is systematically lower in cluster with higher ve- 



locity dispersion. This result would imply that T increases 
with cluster mass, if the shape of the luminosity function 
of cluster galaxies were independent of cluster mass. 

These results are in agreement with the mass-to-light 
function of the ACDM model, but not with that of the 
rCDM model. 

On the largest scales, the behavior of the mass-to-light 
function is less clear. Virialized objects of supercluster 
mass are expected to be extremely rare — there are cer- 
tainly none in the NOG volume. It is possible to obtain 
direct estimates of the global mass-to-light ratio can ob- 
tained from studies of peculiar velocities, which probe de- 
viations from the uniform expansion of the Universe. Note 
that the peculiar velocity of a given galaxy is generated by 
a large number of massive halos over a range of distances, 
in both overdense and underdense regions. Thus the mass- 
to- light ratio determined by such studies is an average over 
a range of mass scales. Hudson (1993, 1994) studied the 
density field and peculiar velocity field in a volume very 
similar to the NOG. He assumed a constant mass-to-light 
ration for all systems and found that the average mass-to- 
light ratio of optical galaxies was 30% of the mass-to-light 
ratio of a critical density Universe, if galaxies are unbi- 
ased tracers of the mass. It would be interesting to see 
how these conclusions would be altered if one adopted the 
mass-to-light function found here for the ACDM model. 

Our results can also be used to determine the efficiency 
of star formation, namely the fraction of all baryonic ma- 
terial which has turned into stars. We can convert our 
T function into a stellar-to-total mass ratio if we assume 
a stellar mass-to-light ratio averaged over spheroids and 
disks, of T st ar = 3.4T Q (Fukugita, Hogan & Peebles, 1998) 
in the B-band. Near the minimum of the T function, this 
yields a mass fraction of stars, M sta r/^tot ~ 0.035/i^ . 
The mass fraction in stars can then be converted into the 
fraction of baryons in stars, if we assume that the baryonic- 
to-total mass fraction in virialized systems is the same as 
the global value, £lbaryon/^m- Using f2 b aryon = 0.022/i^q 
from Netterfield et al. (2001), we obtain the stellar- to- 
baryonic fraction 

M star /M baryon ~ 0.26 (^^) 1 (£f ) ^ (7) 



Thus in those dark matter halos which have luminosities in 
the range 1 — 10L* al , masses in the range 10 12 — 10 13 m Q 
and where T ~ IOO/175Y0, ~ 25% of the mass is con- 
verted into stars. This is considerably more efficient than 
the global stellar fraction of <; 10% (Fukugita et al. 1998). 

For the ACDM model, the mass-to-light-ratio increases 
to T = 350/175X0 at cluster masses ~ 5x 10 14 /i 7 ~ 1 m Q . 
At this mass, the star formation efficiency has dropped 
to 7.4%. If we also allow for the redder population of 
such systems, and assign a higher stellar M/Lb = 4.5 
(Fukugita ct al. 1998), the star formation efficiency would 
be 10%. This is in agreement with the results of Balogh et 
al. (2001). For the rCDM model, the mass-to-light ratio 
is ~ 33O/175Y0 for all systems with masses greater than 
lO 13 /i 7 ~ 5 1 m0, and in this mass range the star formation ef- 
ficiency would be ~ 25%. 
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Fig. 6. — The number of galaxies with luminosity greater than 
L* al hosted in a halo is plotted as a function of the system luminos- 
ity. Bars represent ±1<t errors calculated as the standard deviation 
in each bin of absolute magnitude. The mass scale is obtained using 
the mass-to— light ratio derived in a ACDM cosmogony in §4. Also 
shown the best fitting curve as determined in equation refill. 

Finally, we can compare our VSLF to the predictions of 
semi-analytic models, such as those of Benson et al. (2000). 
In Figure 5, we show our LF and that of MFW in compar- 
ison with the predictions of their ACDM and tCDM mod- 
els. Benson et al. constrained their models to match the 
amplitude of the ESP galaxy LF in the £?j-band (Zucca et 
al. 1997) at their i* al . Their tCDM model is quite a poor 
fit, failing to match the slope of our VSLF at intermedi- 
ate luminosities and predicting too many low-luminosity 
objects. Their ACDM model provides a marginal fit to 
our results: the bright-end slope is correct, although the 
predicted LF slope in the range of groups is somewhat too 
steep. This arises because their T function increases as 
m ' 3 between masses of 10 12 mQ and 10 14 m©, whereas we 
find that the behavior of T is quite flat in this range (see 
Fig. 4). An increase in T from the scales of large galax- 
ies up to clusters is also predicted by the galaxy forma- 
tion models of Kauffmann et al. (1999), with <x m 02 
over the range 10 12 TOq < m < lO 15 m0 for both their 
ACDM and rCDM models. This is similar to that found 
by Benson et al., and suggests that the Kauffmann et al. 
VSLF would be qualitatively similar to that of Benson et 
al. The T functions obtained from the semi-analytic mod- 
els of Somerville et al. (2001) also increase with increasing 
mass, however their T values are larger than those of Ben- 
son et al. and are thus in poorer agreement with those 
derived here. 

4. HALO OCCUPATION NUMBERS 

For many purposes, it is important to know the number 
of galaxies which populate a given halo. Once this function 



is known, one can use the well-understood clustering of ha- 
los to obtain the clustering of galaxies (Neyman, Scott and 
Shane 1953; Kauffmann et al. 1999; Diaferio et al. 1999; 
Benson et al. 2000; Seljak 2000; Peacock and Smith 2000; 
Scoccimarro et al. 2001; Benson 2001). 

In Fig. 6, we plot the number of galaxies with luminos- 
ity greater than L* al as a function of the corrected halo 
luminosity and of the total halo mass derived under the 
assumption of a ACDM cosmogony (§3). Note that NOG 
sample is complete for L ga \ > L* al (see Fig. 6 of Paper 
V) and no luminosity selection effects are polluting the 
observed trend. Note that, by definition, we are only con- 
sidering the high- luminosity (L ga \ > £* al ) end of the T 
function. We find that the number of galaxies above the 
Z/* al luminosity threshold, scales with the total luminosity 
of the system as follows 

(\ 0.83±0.25 
jtH . ( 8 ) 

for M s — 51og/i75 < —21. Using the mass-to-light ratio 
derived in section §3 we can derive the dependence of the 
number of objects residing within halos in terms of the 
halo mass and of the specific cosmological model adopted. 
We find that in a ACDM and tCDM cosmogonies, the 
scaling of the halo occupation number with mass is given 
by the following expressions 

/ \ 0.55±0.26 

iV gal (>L* al ) = 6.3^ 3 .10- 8 ^j (9) 
for m ;> 5xlO 12 ft-7 5 1 m0 and 

/ s 0.83±0.33 

N 9 al(> L* g J = 4.2h 75 ■ 10- 12 — (10) 

\m Q j 

for m ^ 10 13 /iy 5 1 m©, respectively. 

Peacock and Smith (2000) performed a similar analysis 
using CfA (Ramella, Pisani & Geller 1997) and ESO Slice 
Project groups (Ramella et al. 1999). Their Figure 6 sug- 
gests power-law exponents of ~ 0.55 and ~ 0.7, for ACDM 
and rCDM models respectively, in good agreement with 
our results. 

5. THE RELATION BETWEEN OPTICAL AND X-RAY 
PROPERTIES IN VIRIALIZED SYSTEMS 

Having derived the functional form of the optical lumi- 
nosity distribution of virialized systems, we now investi- 
gate the relation between optical and X-ray emission prop- 
erties in galaxy systems. The techniques involved are the 
same as those developed in section §2. The general prob- 
lem of deriving the functional form Lx = f(L), relating 
the optical and the X-ray luminosities, can be formulated 
using the following differential equation 

r £x = f(£s) 

I (ii) 

I df(L s ) _ <f(L s ) y ' 

[. dL s 0(L X ) 

where 4>(L S ) is the optical LF of systems emitting in the 
X-ray band with a luminosity distribution (/>x(£x), an( l 
where Cx = /(£ s ) is the initial condition, which can be 
determined from the data. 
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Fig. 7. — Upper: the local X-ray luminosity functions in four 
different bands are compared over a luminosity range that de- 
scribes both groups (faint end) and clusters (bright end). The 
shape of the NOG VSLF is shown for comparison. Lower: the 
ratio f between the logarithms of the four X-ray luminosity func- 
tions and a power law expression of parameters ctx = —1.85 and 
0*. = 4 • 10- 7 Mpc- 3 (10 44 crgs- 1 )«x-l i s shown as a function of 
luminosity. 

We assume that all virialized systems detected in the op- 
tical also radiate in X-rays. Their optical LF is a Schechtcr 
function whose parameters (<fi*,a s ,M*) are known over 
the luminosity range —18 < M s — 51og/i7 5 < —25 (§3 and 
Fig. 7). The XLFs in various bands are given by Ebeling 
et al. (1997). 

The optical emission is expected to be connected with 
the X-ray emission under the standard assumption that 
the same gravitational potential shapes the gas density 
distribution as well as the galaxy distribution (Cavaliere 
& Fusco-Femiano 1976). We can gain some insights into 
the optical-to- X-ray luminosity ratio making some approx- 
imations. Fig. 7 shows that (f>x(Lx), for both groups and 
clusters, can be well approximated, in almost all the X- 
ray bands, by a single power-law expression 4>x{Lx) = 
4>x L ax . In particular, in the bolometric and 0.1-2.4 keV 
bands, the local XLFs of Ebeling et al. (1997) are very 
well approximated (and extrapolated into the group do- 
main) by the following fitting parameters (fr x ~ 4 • 1CP 7 
Mpc" 3 (10 44 erg s _1 ) ax_1 and a x ~ 1.85. In this case, 
the particular solution of eq. 1 1 is simple and can be ex- 
pressed by the following analytical formula 



L X = {f(£ s f* C(3x [l\ (A,, ±f) - T 4 (ft, £)] } /<x 

(12) 

where C — </>* (3 — 1 + a (for s and X subscripts) and 
Fj is the incomplete gamma function. However, if we in- 
tegrate numerically the differential equation 11, we obtain 
the result plotted in Fig. 8. 




Fig. 8. — The relation between optical (L s ) and X-ray (Lx) lumi- 
nosities as determined in eq. 12 is shown together with the observed 
luminosities of Abell clusters. The straight line is the scaling rela- 
tion predicted in a ACDM universe using the optical (this paper) 
and the X-ray (Ledlow et al. 1999) mass— to-light ratios with arbi- 
trary normalization. The thick solid line is the exact solution and 
the dotted line is the analytical solution of cq. 12. The hatched 
region corresponds to the errors in the Lx — L s relation obtained 
from the errors in the respective luminosity functions. 

In order to test the solution and determine the normal- 
ization constant, we have used a sample of Abell clusters 
for which both the optical blue magnitudes and the X- 
ray luminosities in the 0.1-2.4 keV band are known. The 
largest local samples of X-ray luminosities compiled to 
date are the X-ray Brightest Abell Cluster Sample XBACS 
(Ebeling et al. 1996) and the Brightest Cluster Sample 
(BCS) (Ebeling et al. 1998). We found that 28 XBACS 
and 3 BCS clusters also have estimated Bj luminosities 
(Girardi et al. 2000). At the low- luminosity end we used 
a sample of compact groups whose luminosities in the X- 
ray (bolometric) and optical (blue) passbands are given 
by Ponman et al. (1996) and Hickson et al. (1992) respec- 
tively. In order to reduce observational errors, the initial 
condition has been determined by averaging the luminosi- 
ties of the four optically brightest clusters in the sample. 

The overall agreement between our prediction and the 
data is shown in Fig. 8. Note that the Girardi et al. (2000) 
data shown here are at the tail of our optical VSLF. Of 
particular interest is the apparent break in the power-law 



-21 +5 log /i 75 ,Z, x ~ 10 42 erg s^h 



75 



behavior at M s 
Fainter than this break, the numerical solution is well ap- 
proximated by Lx oc ig 5 , whereas in the brighter regime, 
both numerical and analytic solutions are well approxi- 
mated by Lx oc Lf 5 . 
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When considering X-ray properties, as well as mass-to- 
light ratios, it is convenient to separate poor systems, by 
which we mean systems with one or a few L* al galaxies 
(poor groups), from rich systems (rich groups and clus- 
ters). We define the former to be those systems with 
-21 > M s - 5 log h 75 > -23 or L x < 10 42 erg s^h^ 2 , 
and the latter to be systems with —23 > M s — 5 log /175 > 
-25 or 10 42 < Lx/(erg s" 1 /^ 2 ) < 5 • 10 44 . Beyond this 
range the X-ray and optical LFs are not well determined 
from our data 

6.1. Scaling relations for rich systems 

In the regime of rich systems, as defined above, the X- 
ray-to-optical scaling is Lx oc L 35 . In this regime, for a 
ACDM cosmology, we found in §4 that 

to oc Ll- 5±0 - 26 . (13) 

This then implies that to oc l^ 43±0 08 . This is in good 
agreement with the analysis of Ledlow et al. (1999), who 
compared the observed XLF to PS mass functions and 
found to oc i^ 4±0 03 , in the bolometric X-ray band over 
the range 10 41 < L x /(erg s" 1 /^ 2 ) < 5 • 10 45 ) (solution 
ACDM2 in their Table 1 which is similar to the ACDM 
model adopted here.) 

It is possible to compare this scaling with other cluster 
parameters such as velocity dispersion a and temperature 
T. Under conditions of spherically-symmetric and isother- 
mal equilibrium, the X-ray luminosity is connected to the 
velocity dispersion of the virialized halos by the relation 
Lx oc f 2 cr 3 T^ 2 (Quintana & Melnick 1982), where / is the 
ratio of the gas mass to the total cluster mass. If the gas 
and galaxies are in hydrostatic equilibrium, the average 
plasma temperature is proportional to the depth of the 
cluster potential well (T oc cr 2 ), leading to the following 
simple scaling relation Lx oc f 2 cr 4 . The observed Lx — a 
relation for clusters is not very far from this theoretical 
prediction. Several authors reported an empirical Lx — o 
relation close to Lx oc cr A (Quintana & Melnick 1982; 
Mulchaey and Zabludoff 1998), while Xue & Wu (2000) 
and White, Jones & Forman (1997) derived somewhat 
steeper relations, L x oc <j 5 - 3±0 - 21 and L x oc cr6- 3 8±o.46_ 
respectively. 

Our results allow an alternative determination of this 
scaling relationship. Since our systems are defined by 
a fixed overdensity criterion, their typical size scales as 
R s oc [J*^° L^(L)dL}~ 1/3 . Using eq. 5 of paper V and 
the cosmology dependent relation to = L a we obtain, 

L s oc 0-3^1 . (14) 

Using the observed scaling relationship for the X-ray-to- 
optical luminosity ratio (Lx oc L% ) we can write 

60 

L x oca-3°-i (15) 

For rich systems in the ACDM cosmology, this yields 
L s oc fj 1 - 7 ^ - 38 which is in good agreement with the re- 
sults of Schaeffer et al. (1993), who found L v oc f7 1 -87±o.44 
and those of Adami et al. 1998 who reported Lb } oc cr 1,56 
with a large scatter. 

In the X-ray band, our result would translate into the 
following scaling relation, Lx = cr e±1 - 3 , consistent with 
those found by Xue & Wu (2000) and White, Jones, & For- 
man (1997). In rCDM cosmology we would have obtained 



the following relations L s oc cr 3±0 13 and Lx = a 10 - 5±1 - 3 , 
in disagreement with the observations. 

We can go a step further if we assume that clusters fol- 
low an isothermal relation a oc T 1 / 2 (which appears to 
be supported by the data of Ponman et al. 1996; White, 
Jones, & Forman 1997; Xue & Wu 2000). Inserting this 
a — T model into the Lx — & scaling relationship we ob- 
tain Lx oc T 3±0 - e5 . Again this result is in agreement with 
the fit of White, Jones, & Forman (1997) and Xue & Wu 
(2000). 

Thus, in the regime of rich systems, the derived scal- 
ings are in agreement with the observations for the ACDM 
cosmology. For the tCDM cosmology the agreement is 
poorer. 

6.2. Scaling relations for poor systems 

For poor systems, i.e. super-L* al galaxies and poor 
groups, we found Lx oc L 4 ' 5 . In this regime the mass-to- 
light ratio is nearly constant, m oc L 4±01 , for both ACDM 
and rCDM cosmologies. This yields L x oc L 2 2±0 - 3 , which, 
besides being in agreement with what we found, confirms 
in an independent way the bivariate scaling behavior of the 
optical-to-X-ray luminosity ratio over the group and clus- 
ter scales. Extending this analysis to velocity dispersions 
using 14, we derive a steep relationship L s oc cr 3±0A5 . Now 
coupling our results with the observationally determined 
Lx — L s scaling relationship we obtain Lx oc (j 3 - 9 ^ 1 - 2 . 

Because of the problems mentioned in §1, direct obser- 
vational results for groups and poor clusters are much scat- 
tered, although they tend to lead to a flatter Lx — cr rela- 
tion compared to that of clusters. Mulchaey & Zabludoff 
(1998) obtain Lx oc cr 4 - 3 ±°- 4 7 p on man et al. (1996) ob- 
tain Lx oc cr 4 9±21 , Helsdon & Ponman (2000) give Lx oc 
o- 4 - 5±11 , whilst Mahdavi et al. (1997) and Mahdavi et al. 
(2000) found much flatter relations, i.e. Lx oc a 1 - 5 ^ - 25 
and Lx oc (jO.4io.3_ Thc fit of Xue & Wu ( 2000 ) ig intcrmc _ 

diate in slope with Lx oc cr 2 - 35 ± - 21 . Our derived relation 
is consistent with the steep slopes and is marginally con- 
sistent with that of Xue & Wu. 

Now if we assume a oc T (Helsdon & Ponman 2000), 
we obtain the steeper relation Lx oc T 3 - 9±1 - 2 which is in 
agreement with the theoretical predictions of Cavaliere, 
Menci & Tozzi (1997, 1999) (see the discussion at the end 
of the section) and in excellent agreement with the results 
Lx oc T 4 9±0 - 9 of Helsdon & Ponman (2000). 

If we had assumed a condition of virial equilibrium 
between condensed and diffuse baryons on the scale of 
groups, we would have obtained a flatter relation Lx oc 
T 5 similar to the trend observed in clusters. Since this 
does not provide a good fit to data (see Fig. 9), we rein- 
force the conclusion of Ponman et al. that at the low end 
hierarchy of the clustering pattern, i.e. for low X-ray tem- 
perature systems, the virial equilibrium condition between 
galaxies and gas does not apply. 

6.3. Discussion 

The break in power-law scalings of mass, Lx, cr and 
T between rich and poor systems, occurring at Lx ~ 
10 42 ergs l hj§, has been noted by many authors. Here 
we have found that a change in the shape of the Lx — L s 
relation occurs at the same point. 
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Fig. 9. — The relation between the bolomctric luminosity and 
the temperature in ACDM cosmology. Data are from Arnaud & 
Evrard (1999) (filled squares), Allen & Fabian (1998) (triangles) 
and Ponman et al. (1996) (empty squares). The thick line refers to 
our predicted scaling behavior. The dashed line refers to the self- 
similar case and the solid line to predictions of the thermodynamic 
model of the z = ICM (in a ACDM scenario) developed by Tozzi 
& Norman (1999) (the entropy excess is K = 0.3 ■ 10 34 ergcm 2 g~ 5 /» 
constant with epoch and the cooling is included). The dotted line 
shows the relation Lx — T at z=0 derived by Tozzi and Norman 
(1999) within the projected radius used in Ponman et al. (1996). 

At a similar point there is a net change in the chemical 
properties and the spatial distribution of the ICM on the 
scales of groups (Rcnzini 1997, 1999). One might think 
that this is due solely to a change in the properties of the 
ICM. However, if the ACDM PS mass function is correct, 
then it is interesting that this transition corresponds to the 
same mass scale (M s ~ —23 + 5 log h 75) which separates 
efficient galaxy formation, with T low and approximately 
constant, from inefficient galaxy formation at which T be- 
gins to rise as I/°' 5 . This would suggest a connection be- 
tween galaxy formation and the properties of the ICM. 
There are two proposed models for this connection. The 
star formation may have preheated the gas, or the low- 
entropy gas may have preferentially cooled into stars. 

The broken power-law behavior is in quantitative agree- 
ment with predictions of the preheating scenario (Kaiser 
1991, Evrard & Henry 1991) in which the entropy of the 
hot, diffuse intracluster medium is raised at early time 
(prior to gravitational collapse) by non-gravitational heat- 
ing such as feedback effects of star formation, SN winds, 
shocks etc. Cavaliere, Menci, & Tozzi (1997, 1999) pro- 
posed a semi-analytic model of shocked ICM gas to ex- 
plain the the observed Lx — T relation. The same re- 
lation is reproduced by a model of Balogh et al. (1999) 
using the physics of an adiabatic isentropic collapse of 
a preheated gas. Both models make detailed predic- 
tions for the Lx — T relations which are approximated 
by the power laws Lx oc T 5 in the luminosity range 
Lx < 0.25 • 10 44 (erg s" 1 /^ 2 ) (groups) and L x oc T 3 in the 



luminosity range 0.25 • 10 44 < L x /(erg s" 1 /^ 2 ) < 5 • 10 44 
(rich groups and clusters). Tozzi & Norman (2000) com- 
bine these two scenarios in a single model taking into ac- 
count an initial entropy excess and the transition between 
the adiabatic and shock regime in the growth of X-ray ha- 
los. Their Lx — T model and our predicted scalings are 
compared to data in Fig. 9. 

On the other hand, the preferential cooling of low- 
entropy gas into stars in galaxies might also break the 
scaling (Thomas & Couchman 1992; Bryan 2000; Muan- 
wong et al. 2001). For this process to explain the scaling, 
one requires greater efficiency of star formation for low- 
mass halos, e.g. 10% at T ~ 1 keV vs 4% at T - 10 keV 
(Bryan 2000). As discussed in §4, for the ACDM mod- 
els, we found that the maximum star formation efficiency 
was as high as 25% for poor systems, about 2.5 times the 
global efficiency of 10%. For rich systems the star for- 
mation efficiency drops below the global average. Note 
that for a different mass spectrum, e.g. tCDM there is no 
mass dependence of the star formation efficiency for high 
masses. 

In summary, at least for the ACDM spectrum, it is in- 
teresting that the break in X-ray scaling properties occurs 
at roughly the same mass as the change in star formation 
efficiency. This suggests that there is a physical link be- 
tween galaxy formation and the X-ray properties of the 
hot gas. Of course increased galaxy formation will lead to 
greater feedback as well as more efficient cooling of low- 
entropy gas, so it is difficult to distinguish between the 
scenarios discussed above. 

7. SUMMARY AND CONCLUSIONS 

We have used the B-band luminosity function of viri- 
alized systems to investigate scaling relations over a wide 
dynamic range of mass, from single galaxies to clusters of 
galaxies. 

When our LF is compared to the Press-Schechter mass 
functions predicted in CDM cosmogonies we find that all 
these models fail, if a constant mass-to-light ratio is as- 
sumed. Specifically, in the tCDM and ACDM models, 
the mass-to-light ratio must vary as L -0 5 to match the 
faint-end of the luminosity function. On the other hand, 
a constant mass-to-light ratio and a ratio varying as L 5 
match the bright ends of the rCDM and ACDM models, 
respectively. For the ACDM model, the efficiency of star 
formation is maximized in the range of one to ~ 5 L* al 

or lO 12 5 to lO 13 5 /!^ 1 m©, the regime of single galaxies 
and poor groups. In this range, the fraction of baryons is 
stars is ~ 25%. The latter behavior of the mass-to-light 
ratio is in qualitative agreement with the predictions of re- 
cent semi-analytical models of galaxy formation, in which 
galaxy formation is inhibited by the reheating of cool gas 
on small-mass scales and by the long cooling times of hot 
gas on large-mass scales. The variation of the mass-to- 
light ratio with scale found by some authors through direct 
estimates of the masses of galaxy systems lends indirect 
support to the ACDM cosmology. A further test of the 
varying mass-to-light ratio model could be made via the 
comparison of observed peculiar velocities with the pre- 
dictions from nearly all-sky optical galaxy catalogs such 
as the NOG. 

Since our sample of galaxies is complete for objects 
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brighter than L* al , we have also measured the halo oc- 
cupation number. We find that, for ACDM models, this 
quantity grows more slowly than the total mass of the sys- 
tem, as required in the "halo" model of explain clustering. 

Comparing X-ray luminosities as a function of optical 
luminosities in virialized systems, we find a break in the 
power-law slope of the relation as we go from groups to 
clusters, independent of the cosmological model. Dynam- 
ical scaling relations in X-ray systems are known to have 
break at a similar mass scale. Furthermore, the break oc- 
curs at a similar mass scale to that at which the mass- 
to-light ratio changes in the ACDM model. This sug- 
gests physical link between galaxy formation and the X-ray 
properties. However, these data cannot say whether this 
is due to reheating of the hot halo gas by supernovae or 
efficient cooling of low-entropy gas. 

We look forward to forthcoming large field redshift sur- 



veys planned with SDSS (Gunn & Knapp 1993), 2dF (Col- 
less 1998), 2MASS (Huchra et al. 1998), which will (a) 
extend the range over which the system luminosity is de- 
termined (b) lower the luminosity threshold over which 
group members are selected and (c) reduce the errors in 
their distribution statistics. Moreover the DEEP2 redshift 
survey (Davis & Faber 1998) will probe the variation with 
cosmic time of the M/L — L, Lx — L and L — a relations, 
allowing a deeper understanding of galaxy formation. 
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Masses of systems for different models 
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